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ABSTI^ACT 


This work deals with the themriomechanical response of functionally graded 
matenals A plate like structure subjected to sharp temperature gradient m 
the thickness direction is considered Theoretical models and numerical 
results for temperature field, heat flux, thermal residual stress, and thermal 
stress are presented both for functionally graded and equivalent bi-layered 
materials Significant trends and design guidelines to select appropnate 
composition profile for specific material pairs are presented. The role of 
microstructure, especially m regard to percolation phenomenon m transport 
properties and its significance for functionally graded matenals is discussed 
and a novel approach to incorporate such effects m design calculations is 
presented. 
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CHAPTER 1. 


Introduction 


1.1 Need for Functional gradient material (FGM). 

In recent years the environments in which matenals are used have 
become more demanding. Conventional materials cannot withstand the 
severe environments confronting modem technologies. Further 
developments in science and technology rely heavily upon the 
development of new materials that can withstand the conditions that are 
created by advancing technology 

In the development of these new materials there are two approaches. 
One IS to develop new material that differs completely from any present 
matenals The other option is to develop new function for existing 
matenals The development of composites satisfies both of the above 
approaches. Many new types of composite have already been fabncated 
in accordance with the material properties demanded by today’s 
teclmology. 

A general requirement for industrial materials is uniformity m 
matenal properties That is, for mdustnal applications, it is essential that 
every part of the matenal in use exhibit uniform properties. The 
requirements for conventional composites aie no exception. Most effort 
in the development of eomposites in the past has been put into 
determining how to uniformly mix the dispersoid within the matrix. As 
a result, from macroscopic viewpoint, even composites were regarded as 
homogeneous because in these matenals the mechanical properties and 
other material characten sties prove to be homogeneous 

In contrast, studies are also being conducted to design materials that 
have two different functions within the given matenal These materials 
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are inhomogeneous composites that are charactenzed by having 
different material characteristics on separate surfaces or m separate 
parts An example of a composite having different functions in its 
different parts is a coated or joined material designed to improve a 
material’s suifacc chaiactciistics However, these inhomogeneous 
composites possess sharp boundaries The boundary often exhibits 
vanous undesirable pioperties caused by the existence of discontinuities 
in the mateiial’s mechanical, physical, and chemical characteristics at 
the boundary An example is a separation at the boundary due to thermal 
stress 

For this reason a proposal to design a new material aimed at 
eliminating the macroscopic boundary m lammatcd-type material is 
being examined 'fliat is, the appioach is to synthesize inhomogeneous 
composites in which, the material’s mechanical, physical and chemical 
property changes continuously, and which have no discontinuities 
within the material These materials are called Functional Gradient 
Matenals (An alternative name is FGM). Fig 1.1 shows schematically 
the concept of a Functionally Giaded Material 

1.2 Definition of FGM and their importance. 

FGM are those materials, which exhibit a progressive change in 
composition, structure, and properties as a function of position within 
the material See Fig 1 2 

FGM are designed to take advantage of the attractive features of each 
of Its constituents For example, one constituent may be a ceramic, 
which offers good high temperature behavior but is mechanically bnttle. 
Another may be a metal which exhibits better mechanical and heat 
transfer properties but cannot withstand exposure to high temperature A 
structural component exposed to regions of both high and low 
temperature could thus be predominantly ceramic within the hotter 
region and metal within the cooler region. Sec Fig 1 .3. 

The design of an b’GM is thus centered on the optimization of its 
structure in order to maximize (or minimize) some property associated 
with Its in-service performance, while simultaneously maintaining other 
properties within acceptable limits. 

Another aspect lies m predicting the charactenstics of an FGM, for a 
given composition profile, during fabrication and under in service 
condition Use of theoretical models to aid in FGM design and to predict 
FGM fabncation and m-service behavior is of crucial importance. 
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I'he fact that the composition of an FGM can vary over such a wide 
range means that a variety of fundamentally different microstructures 
can exist across the graded direction This, m turn, means that the 
thermo-physical properties, which are generally strongly dependent on 
the microstructuie, will also vary with position within the material A 
realistic model must appiopi lately account for this fact 

1.3 Approaches to obtain FGM. 

These are two approaches to obtaining FGMs, as shown m Fig 1.1. The 
first method is to eliminate the boundary of laminated-type composites 
[from (c) to (b)j, thcicby eliminating discontinuities in the properties at 
the boundary 'fhe second method is to make nonuniform distribution of 
dispersoids in homogeneous composites [fiom (a) to (b).], thus creating 
multiple functions within the matciial 'Ihe giadicnt can be assigned not 
only m the diiection of material’s thickness but in the direction of its 
width as well 

In nonnal FGMs, the concentration of the dispersoid changes 
gradually within the composite, but another type of FGM can be 
designed by continuously changing the morphology of the dispersoid, 
state, crystal structure and distnbution pattern (Hiran, 1996). 

1.3.1 Continuously changing the morphology of the dispersoids. 

Spherical and fiber like dispersoids are distributed over surface X and 
Y. The shapes of these dispersoids are gradually varied (Fig 1 3a), 
Within a matenal, the dispersoid will be changed from spherical to 
fiber-like In this way a high Young’s modulus and a high toughness 
value can be obtained on surface X while also achieving a high strength 
on surface Y, If surface Y is changed to a composite with a flake-like 
dispersoid it will be a material that is rich in lubricity. Plant systems 
observed in trees and bamboo often have this type of structure(MRS 
Bulletin, 1995) 

1.3.2 Continuously changing the state. 

With dense materials on surface X, introducing small voids into that 
surface can create a porous surface Y. Bioactive materials are an 
example of this type of composite The advantages are that surface X 
has sufficient strength to be used as a structural material, while surface 
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Y has good compatibility necessary for a biomaterial These voids also 
exhibit superior thermal stress relaxation This characteristic can be used 
to design unique and intriguing FGMs (Hirari, 1996) 

When controlling the level of sintering from the front surface to back 
surface continuously changes the void distribution, the result is a 
simultaneous continuous change in the Young’s modulus from the front 
to back This phenomenon is used to assign gradients m the 
piezoelectnc properties within the matenal in order to generate a 
flexural vibration mode (Kawasaki and Watanabe, 1990) 

This technique of continuously changing the state [Fig 1 3b] has some 
interesting applications m polymeric materials For example, by 
gradually lowering the degree of polymerization of a rod-like polymenc 
matenal, starting at the outer perimeter and continuing toward the 
center, it is possible to control the damping characteristics or vary the 
optical properties of the final matenal (Hirari, 1996). 

1.3.3 Continuously changing the crystal structure. 

In as amorphous material, high mechanical strength can be obtained 
by crystallization If a surface X is crystallized while the opposite 
surface, Y, is maintained in the amorphous state, surface X will have 
high mechanical strength while suiface Y will have a good corrosion 
resistance, a gradient of these two properties will develop between the 
two sui faces [Fig 1 3c|. I'oi melallic malciials, an impiovcnieiit in the 
surface properties is obtained by introducing ions to the surface, thereby 
making it amorphous These are examples of crystalhne/amorphous 
FGMs If a crystallographic axis on surface X is oriented m one 
direction, and no specific orientation is given for surface Y, then it is 
possible to effectively utilize these anisotropic properties in order to 
influence the matenal’s properties (Hirari, 1996). 

By continuously changing the crystal structure m the thickness 
direction of a film from f c.c (surface X) to h.c.p (surface Y), an attempt 
has been made to prepare a high density perpendicular magnetic 
recording medium (Osaka et al , 1 990). 

By observing the change from paramagnetic to ferromagnetic that 
results from the martensite transformation, a magnetic FGM that has a 
gradual change m its saturation magnetization can be designed. This is 
achieved by continuously changing the transformation level from 
surface X to surface Y using the matenal to control deformation 
(Watanabe et al , 1993). 
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1.3.4 Continuously changing the distribution pattern. 

High strength can be obtained by uniformly distnbutmg the 
dispersoids on surface X High electrical conductivity can be obtained 
by distributing the dispersoids in the form of a network pattern on 
surface Y [Fig 1 3d] 

1.3.5 Continuously changing the grain boundary characteristics. 

By maintaining a nearly perfect gram boundary on surface X while 
introducing voids or faults into the boundary regions on surface Y, it is 
possible to control the diffusion of the atoms or molecules and thus 
obtain specialized charactenstics suitable for a variety of sensor 
application, or for use as a catalyst 

1.4. Applications of FGM. 

Some important applications of FGMs are given below. 

1. Joining Media 

Since they have a potential for reducing thermal stress, FGMs have 
been investigated for joining metals to ceramic (Metal/FGM/Ceramic) 
and both metals and ceramic to FGMs (Metal/FGM, ceramic/FGM) 

An FGM filler of Ni-MgO was used to join Ni to MgO, forming 
Ni/Ni-MgO(FGM)/MgO. The fracture strength (three-point flexural test) 
of the joints formed using the FGM Ni-NiO (72-128 Mpa) and FGM Ni- 
MgO (110 Mpa) was higher than when Ni was directly joined to MgO 
(30-60 Mpa). Iron was also successfully joined to AIN via an FGM of 
Al-AIN by containing the aluminum surface of the FGM with the iron 
and heating at 650”C for 48 hours in flowing argon. The micro-hardness 
of the Al-rich zone in the FGM was greater than pure aluminum. 
[Atarashiya, 1992] 

2. Shuttle Thermal Insulation. 

A shuttle protection FGM type tile developed for the space shuttle is 
more efficient and less costly than tiles currently being used. This tile, is 
known as Toughened Um-Piece Fibrous Insulation (TUFI), is a low- 
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density composite thermal insulation It is the first FGM composite 
where the density of the material varies from high at the outer surfaces 
to low in the intenor insulation 

The current tiles are a rigid glass-fiber composite and are about 93% 
air, with a thin glass coating on top The Reaction-Cured Glass (RCG) 
coating IS physically much like window glass and is only about 0 03 mm 
thick. Because of the coating gets little support from the underlying tile. 
It cracks or chips easily Unlike RCG, TUFI permeates the pores nearer 
the surface of the insulation material This supports and reinforces the 
outer surface, which makes the surface material less subject to impact 
damage (Aerospace Engr , 1 994) 

3. Surface Wave Devices. 

Tam and Liu (1993) studied a piezoelectric FGM plate having 
continuously changing material properties in the direction of the plate 
thickness by wave propagation analysis The wave frequency spectrum, 
the energy propagation speed, the electromechanical coupling constants, 
and the mode shapes were determined by these analyses The results 
showed that the gradient properties can significantly improve the 
effectiveness of surface wave devices 

4. Optical Fibers. 

Kamaswami and Nishizawa (1965) proposed an optical fiber with a 
continuously changing refractive index in order to increase its 
transmission capacity. It has a higher refractive index at the core, while 
its sheath has a lower refractive index This type of optical fiber is called 
a graded index-type optical fiber (Gl-type). The Gl-type fiber has been 
shown to be capable of carrying more than 10 times the information 
compared to the step index-type optical fiber. Historically, Gl-type 
fibers are the first example of a successful application of the graded 
property concept 

5. Implant materials. 

In order that man-made materials can coexist inside the human body, 
it is highly desirable on organically joining the matenal, that the 
boundary region between the human body and the matenal has a 
functional graded structure. 
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The Ti-Al-V alloy is known to be very compatible with the human body 
and IS widely used foi artificial bones and joints, nevertheless bonding 
them to human bone is still a pioblcm However, a new type of implant 
material composed of a I’l-Al-V-typc substrate has been made. This 
material is coaled with Ti-Al-V beads to make it porous The material is 
further coated with the bioactive hydroxy-apatite (HAP) by plasma 
spraying Easier infiltration of the bone tissue thus occurs with time, 
creating a graded structure between the bone and implant material. The 
result is a better bonding between the human body and the inorganic 
material [Oonishi et al , 1994] 

The bending stiength of HAIVTi is 150 Mpa, which is very close to 
that of a human bone This result suggests that the FGM is effective as 
an implant matenal [Watari, 1994]. 

6. Fire Retardant Material. 

A phosphoric acid and amide mixture can be impregnated into wood 
by a high pressure technique When the phosphorus content is graded 
from the surface to the inner structure by phosphorylation, it becomes an 
excellent fire retardant material [Ishihara et al., 1992]. 

7. Building Material. 

Takahashi et al.( 1993) has successfully developed FGM with 
moisture absorption and release functions made by step wise or 
continuously changing the composition of two components-a zeolite or 
ligneous material-based humidity conditioning matenal, capable of 
absorbing and releasing moisture, and calcium silicate-based concrete. 
This FGM can absorb a maximum of approximately 2 5 times more 
moisture, and is particularly good in its initial response to humidity 
increase or decrease. 

In addition this new material hardly changes its dimensions by 
expansion or contraction due to humidity, and it is strong, fire proof, 
frost damage resistant, and quite amenable to shape forming from flat 
board to curved comer material. As a building matenal, this material 
with a gradient moisture conditioning property from its outer wall to 
inner wall, could eliminate construction steps for faster completion of a 
project. 
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8. Future Shuttle and Aerospace Planes 

Several study programs are underway to establish material design 
technology for creating complex-shaped FGMs and introducing FGM 
designs into the fuselage and engines of a reusable shuttle plane and 
other supersonic planes 

In the development of large FGM test samples, a number of basic 
shapes have been adopted to promote the use of FGMs in the space 
plane fuselage and engines These shapes include a semi sphencal shape 
that IS required to create FGMs for use in the planes nose cone section, a 
hollow cylinder to make FGMs for application m the front edges of the 
wings and engines, and a squre plate to manufacture FGMs for use m 
other sections (MRS Bulletin, 1995). 

9. Magnetic Recording Medium. 

Osaka ct al(1990) picpaicd a CoNiRcP f'GM film by electroless 
plating In this film the structure changes from a randomly onented f.c.c 
structure to an onented hep structure m the thickness direction. They 
reported that by changing the composition and the structure of the film, 
a high density magnetic recording medium can be obtained A magnetic 
FGM was prepaied by applying a gradual inhomogeneous deformation 
to SS304 austenite stainless steel (Fe-18Cr-8Ni) A magnetic FGM can 
be used as a position measuring device by combining it with a magnetic 
sensor For examples, it can be used as a device that determines the 
focal point in an automatic focussing camera [Watanabe et al., 1993]. 

1.5. Present Work. 

The preceding sections have given a bnef overview of the vanous 
form of FGM and a range of current and potential application of FGM. 

The present work is concerned with FGM comprising of two materials 
with a composition gradient along a specified direction in a plate-hke 
structure. Typically, such structures are used for intense thermal 
shielding and arc subjected to sharp temperature gradients. Accordingly, 
metal/ceramic composites with a composition gradient are pnme 
candidates for such applications. Their mechanical integrity under such 
service conditions is of great interest and the emphasis is on specifying 
concentration profile for a given material pair to minimize thermal 
stresses set up in the composite Current approaches to prediction of 
temperature distnbution, thermal stress, thermal residual stress and 
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structure optimization, generally make use of FEM techniques and, are 
reviewed m chapter 2 

The current work describes a simpler analytical/numencal approach 
for determination of tcmpciatuic field and thermal sticss in a FGM and 
these have been described in chapter 3 and chapter 4 respectively 

In chapter 3, major emphasis is given to derive an expression for 
temperature distribution and heat flow rate in FGM, subjected to 
temperature difference across its thickness. Some attention is also paid 
to determine heat flow rate in bimetallic plate, having the same volume 
of each constituent as in FGM Finally, heat flow rate in FGM are 
compared with Bimetallic plate, having same volume of each 
constituent as in I'GM, subjected to same temperatuic diffeicnce across 
their width 

In chapter 4, an expression is derived for determination of thermal 
stress, which develops duiing in-service condition of FGM, due to 
temperature difference across the thickness of FGM Further, thermal 
residual stresses were determined m FGM, subjected to a temperature 
drop of 100°C from a stress free state Thermal stresses are also 
calculated for bimetallic plate, having equal volume of each constituent 
as in FGM 

Determination of temperature field and thermal stresses requires 
models for variation of various physical and mechanical properties with 
composition and suitable mles of mixture have been assumed as per 
general practice. Flowever, the microstructuie of a composite matenal 
may have an important bearing on its transport properties. Accordingly, 
a short discussion on the fabrication methods is also included in chapter 
2 to better appreciate the effect of microstructure on the 
thermomechanical response of FGM. The related issues are discussed 
and a novel approach for incorporating such effects is proposed in 
chapter 5. Major conclusions and suggestions for future work are 
presented in chapter 6 
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(a) - (b) (c) 

f 'y I I Conventional composites and a func- 
Donal gradient malenaL (a) homogeneous composilc, 
(b) functional gradient material, and (c) coated or 
jomed-type composite, . ■' . 



1 ig ! 2 Schciii.ilic illnsliation o (he mitiosliiictuie ol a two-cdmponcnl I'tiM, loi which 
the composition giadually vanes fiom cntiiely the one component to cntncly the other 
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CHAPTER 2 


Literature Review. 

2.1 Methods of Manufacturing Functionally Gradient 
Materials 

The most important objective m the preparation of an FGM is to achieve a 
well-controlled (according to the design) distnbution of composition, 
texture, structure, and other necessary elements. That is, the preparation 
method for an FGM must be able to accommodate an arbitrary number n of 
the compositional distnbution function On continuously changing the 
composition or structural elements within the matenal, the change can occur 
m two directions, m the thickness and width directions Most of the recent 
studies on FGM synthesis have been concentrated on change in thickness 
direction Therefore, m this section the methods used m the preparation of an 
FGM that has continuous change in the thickness direction will be discussed. 
The discussion also includes the method to produce an FGM having a step- 
wise change, as well as a continuous change Laminated-type FGMs, which 
are prepared by first producing sheets of film and then layering the films, 
will also be discussed. 

There are several methods for the preparation of an FGM. However, for 
convenience they will be categorized according to the source matenaTs state 
dunng fabrication of the FGM. They are: gas phase, liquid phase, and solid 
phase approaches 

2.1.1 Vapor Phase Methods 

(a) Chemical Vapor Deposition (CVD) 

The CVD method yields a matenal deposit with source gases as the feed 
stock. This deposit is caused by the application of vanous form of energy 
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(heat, light, plasma, etc ) onto the source gases after they are introduced into 
a reactor. Hydnde, bromide and chloride are generally used for source 
materials When the source material is in liquid or solid form it is vaponzed. 
In recent years, organometallic compounds have often been used as source 
materials By continuously changing the mixture ratio of source gases or by 
controlling the CVD conditions, such as deposition temperature or gas 
pressures the CVD method permits relatively east syntheses of various type 
ofFGMs 

SiC/C FGM coatings were prepared by changing the ratio of S 1 CI 4 /CH 4 
stepwise or continuously m a temperature range of 1300-1500*^0 and a 
pressure range of 1 3-13 kPa on graphite plate, carbon fibers, and a 
cylindrical graphite tube (Ucmuia ct al , 1990). Using the similar conditions, 
an SiC/C FGM of 120 /j/h thickness was coated on a model combustor with 
a cylindrical diameter of 30 mm made of three-dimensional C-C composite 
(Suemitsu ct al , 1993). 

In an effort to impiove the oxidation resistance of C-C composites, C- 
fiber-reinforced SiC sintered bodies, a TiC/SiC FGM layer was coated on 
these substrates. In these cases, in order to minimize the difference between 
the expansion coefficients of the SiC coating and various composites 
matenals, first a layer of TiC is produced on the surface of the composite 
material, followed by a layer of TiC/SiC FGM coating, whose composition 
changed gradually from TiC to SiC This CVD work was performed using 
the T 1 CI 4 -CH 4 -H 2 source gas system at 1350*’C and 8 kPa By changing the 
gas flow rates of T 1 CI 4 and S 1 CI 4 during the deposition piocess, an FGM 
coating of thickness 150 was obtained (Kawai et al., 1990). 

(b) Chemical Vapor Infiltration (CVI) 

The chemical vapor infiltration process, which evolved from the 
chemical vapor deposition for surface coatings, uses porous matenals such 
as ceramic bodies and ceramic cloths which contain many residual pores. 
These materials are first placed inside a CVD furnace as a substrate. By 
decomposition or by chemical reaction of the source gases, a coating is 
deposited on the surface of the open-pore or the space inside of the porous 
substrate. 

Agullo et al. (1993) pre-coated SiC-Nicalon fiber with C/SiC FGM 
(thickness; about 100 mm) by using continuous changes in the gas phase 
composition of tetraethylsilane (SiEt 4 ) and cumene (iPrCcHs) at 770T and 
0.8 kPa. Agullo called it a compositional gradient layer (CGL). This CGL’s 
composition changes continuously from C at the fiber side to SiC at the 
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surface side CVl is used to dcnsify the FGM pie-coated, SiC-Nicalon plain 
weave fabnes above 1 OOO^C 

(c) Physical Vapoi Deposition (PVD) Methods 

In the PVD method, a solid source material is energized using different 
kinds of energy source to obtain vaporized particles, which deposit on a 
substrate to form a film Theie is also a technique to accomplish chemical 
reaction of the vaporized particles in a gas phase The deposition speed of 
the PVD is rather slow Therefore only a thin film can be produced using 
this method However, since this technique does not require high 
temperature heating of the substrate, this method is preferred for the 
preparation of the FGM films used in electronics 

Partially stabilized zirconia (PSZ) with a gradient density of 4.2 -5 g/cm^ 
(70-84% theoretical) fiom the surface to the metal/ceramic interface has 
been prepared on a super-alloy substrate by an electron-beam PVD 
technique (Fntscher and Bunk, 1990) 

2.1.2 Liquid Phase Methods 

(a) Plasma Spraying Methods 

In the plasma spraying method, the spraying source power is transported 
to the plasma jet by a torch nozzle. The molten source material is then 
sprayed on to the substrate to form a coating There are mainly two types of 
plasma sprayed FGM coatings, that is, the porosity-graded coating and the 
composition graded coating Use of plasma spraying in the preparation of 
FGM coatings requires the cleaver design of an apparatus that provides a 
continuously changing mixture ratio of the source power to the torch nozzle. 

The Mg0-Zr02/NiCrAl laminated (7 layers) FGM coating was prepared 
on a steel substrate by a plasma spraying technique using MgO (24 wt.%)- 
Zr02 and NiCr-Al (6 wt.%) powers. The spraying took place in the air using 
a single-torch plasma spray reactor (Eroglu et ah, 1993). 

Fukushima et ah (1990) developed a twin torch for use in FGM 
preparation. In this method two plasma torches are placed so that the 
centerline of each torch is aligned to the point of spray deposition on the 
substrate. Two different source materials are sprayed from each nozzle 
simultaneously onto the same jipot on the substrate to obtain a coating of 
composite. Since each torch can be independently controlled for the required 
spray conditions, this method can easily use two source matenals having 
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widely different melting points, such as metal and ceramic, to form an FGM 
coating Using this technique, an eleven-layered FGM coating of thickness 
0 3-0 4 mm has been prepared using an Ni-based alloy (NiCrAlY) and YSZ. 
By using a smaller spray angle the amount of pores m the coating is 
increased This contributes to improved thermal bamer characteristics. 

(b) Copolymerization Method 

Koike et al (1989) have prepared graded index-type polymer optical 
fiber using the copolymerization method. The source matenals for 
copolymenzation are methylmethaacrylate (MMA) with a refractive index 
of 1.492, and benzyl methacrylate (BMA) with a refractive index of 1.562, 
First a polymer tube is prepared using poly-MMA A monomer mixture, of 
MMA and BMA is then added to this tube and circumferentially heated (60- 
80 ) to obtain polymciization The inner wall of the polymer tube swells, 

interacting with the monomer mixture, and eventually a thin gel phase forms 
on the surface of the wall Then the copolymenzation proceeds towards the 
center of the tube When polymerization the polymerization is completed, a 
polymer solid is obtained that has a radially varied mixture ratio of two 
kinds of polymers from the center to the outer perimeter The resulting solid 
has a radial distribution of refractive indices (Koike, 1991). This graded 
index polymer rod can be heat-drawn at 190-280®C into an optical fiber with 
a graded index 

2.1.3 Solid Phase Methods 

(a) Sintering Methods 

Normal Sintering and High Pressure Smtenng Techniques 

On prepanng FGM by sintering cracks can often appear on the 
sintered body due to differences m the sintering characteristics and the 
mixture ratio of the two source powders Watanbe et al. (1991) controlled 
the shrinkage of the powder mixture by blending fine and course particle 
sizes. 

Rabin and Heaps (1993) prepared six-layered AI 2 O 3 /N 1 FGM using 
the power processing method. Graded compacts were produced by 
sequentially layering the power mixture m the die followed by CIP at 140 
Mpa. I’he sintering was earned out at 1400'’C for 3 h in Ar, or by HP at 
about 1300”C for 1 h at a pressure less than 10 Mpa, or by HIP at about 
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1 350‘’C and 100 Mpa 1 liey concluded that for the best quality FGM, careful 
selection and control of the powder particle size ratio and the preparation of 
a proper compositional distribution profile are necessary 

Temperature Gradient Sintering Technique 

In preparing FGM using source powder having a wide difference in 
their sintering tempeiature (for example metal/ceramic system), one 
sintering temperature cannot be used to obtain a good quality sintered body 
due to the difference m shrinkage characteristics of each source powder To 
overcome this problem, the portion containing more of the higher smtenng 
temperature source is sinteied at a higher temperature and the portion which 
contains more of the lower sintering temperature source must be sintered at a 
low temperature That is, the smtenng must be accomplished under a 
prescribed tempcratuie gradient One way to achieve this temperature 
gradient is to apply additional heat by use of a laser beam or infrared beam 
onto one side of the sample while it is m the sintering furnace. 

It IS also possible to assign a gradient to the density by use of 
temperature gradient sintering Kawasaki and Watanabe (1990) sintered a 
cylindrical PZT power compact in air for 1 h They heated one side of the 
sample by an infrared lamp, and created a temperature gradient of 150°C 
within a depth of 5 mm from the sample surface The resulting sintered body 
was of higher density m the upper portion, while the lower portion was of 
low density. Such a density difference gives a continuous change in the 
piezoelectnc characteristics within the material. 

Plasma Activated Smtenng Technique' 

In plasma activated sintering an instantaneous pulsed electric current 
IS applied to the subject powder to initiate current discharges m the voids 
between the powder particles Successful smtenng can be accomplished by 
use of the heat generated by these discharges. Due to this current discharge, 
the powder surface is purified and activated. This method requires a 
relatively short time for sintering, and thus it is easier to control the gram 
growth. This method is also suitable for smtenng materials with lower 
melting points (Bennett et al., 1968) 

Using this technique a laminated-type (8-layers) FGM of the YSZ/SS410 
system, which has a wide difference in sintering temperature between the 
two sources, was produced In this process the YSZ side was heated to 
1200®C, while the SS side was heated to lOOOV using a specially shaped 
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giaphile resistance healei, in oidei to obtain a temperature gradient (Omori 
et al , 1994) 

(b) Self-Propagating High-Temp eratui e Synthesis (SHS) Methods 

This technique uses an exothermic reaction at temperatures exceeding 
2000-3000 to obtain reaction products at a relatively high speed Because 
of Its high-speed reaction the diffusion of atoms is prevented and thus it is 
possible to obtain a graded composition When pressure (water pressure, gas 
pressure etc ) is applied during the SHS process the synthesis of dense 
composites can be achieved 

Gas Reaction Sintering Technique. 

Ni powder and Al fine powder (0.42 /j/« in diameter) are sintered in 
Ni using the exothermic reaction of AI-N 2 to obtain AIN/Ni FGM. Using a 
similar technique, AIN/Al FGM can also be produced (Atarashiya et al., 
1993). 

Hydrostatic Compression SHS Technique 

T 1 B 2 , Ti, B and Cu powders are used as the source for the preparation of 
T 1 B 2 /CU FGM by this technique (Sata, 1993) Using the automatic powder 
spraying and stacking device, these powders were sprayed onto a Cu 
substrate The stacked compact was then ignited at room temperature under 
a high hydrostatic pressure of 58 Mpa T 1 B 2 /CU FGM with diameter 30 mm 
and thickness 1 mm was obtained on the Cu substrate. 

(c) Martensite Transformation Technique 

Watanabe et al. (1993) have attempted to prepare an FGM using 
crystallographic transformation. The paramagnetic phase in austenite 
stainless steel (Fe-18 Cr-8 Ni) transforms into the ferromagnetic martensite 
phase by plastic deformation The amount of martensite increases with 
increasing deformation. Thus the saturation magnetization of the deformed 
austenite stainless steel increases with increasing deformation (strain). Using 
this phenomenon a magnetic gradient function can be assigned by 
inhomogcncously deforming the stainless steel. 
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2.2 Models for theriiiophysical properties. 


The problem of calculating effective thermophysical properties of 
heterogeneous materials is an old one However with new application for 
such matenals continuing to emerge, interest and activity in this subject 
continues to the present day New approaches are sometimes developed with 
specific applications in mind (e g electrical, magnetic, superconducting, 
mechanical, or thermal properties), and, consequently, are not readily 
accessible to all that might be interested This is unfortunate, because their 
applicability may extend well beyond their immediately intended use Some 
particularly comprehensive reviews, related to this general subject, have 
been published within the past several decades 

Models for specific thermophysical properties of a hetrogeneous 
matenal are based on some kind of microstructure-dependent “mixture rule”. 
Wakishima et al , used a rule of mixture which is discussed here in some 
detail because of its application to the modeling of FGMs. They considered a 
matenal having two components, denoted by A and B. Let Pa and Pb be the 
values of some particular property for pure A and pure B, respectively, and 
let their respective volume fractions be c and c, , where c, =l-c if the 
material is 100% dense. For an FGM, these volume fractions are dependent 
upon position along the giaded direction. The well-known Voigt-type 
estimate for the effective value, P, of this property is 

P = P,c + P,c, (2.1) 


Which IS simply an anthmetic mean On the other hand, the likewise well- 
known Reuss-type estimate is given by 


1 _ c c, 


(2 2 ) 


which is a harmonic mean. As the authors pointed out, these expressions 
have only limited validity. They discussed the use of a more general 
expression 


P = cP^+ c^Pg + cc^Q^ (2.3) 

Where Qab is a function that depends on Pa, Pr, c and C/, as well as on 
microstructurc-dcpcndcnt quantities. 
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A question that must also be addiessed is tliat of how much 
information regarding microstructure must be included m order to describe 
adequately mateiial pioperties and behavioi Becker and Richmond treated 
this particular issue, as it relates to effects of grams, pores, and second-phase 
particles m numerical simulations of the mechanical response of ductile 
matenals. Their models were limited to two dimensions, although the need 
for inclusion of micro-structural features, m order to be physically realistic, 
was demonstrated 

Additional consideration of mixture rules for thermophysical 
properties were given by Hirano et al (1988), who related the pertinent mle 
for various properties to the nature of microstructure The specific 
microstructures they included were laminates, fiber composites with fiber 
aligned transversely and longitudinally to the graded direction, thin-layered 
composites, and composites in which the second phases are either flake-like 
or spherical inclusions. 

Another way to estimate properties of heterogeneous media is to use 
statistical information on the microstructure to compute rigorous bounds on 
property values. Torquato (1991) reviewed these methods and showed that 
these bounds can often be used to provide good estimate of properties 
despite having incomplete statistical information on the microstructure. In 
some cases, these methods can result in improved predictions over bounds 
computed solely from the volume fraction occupied by the constituent 
phases and from individual phase properties 

Still another approach to describing transport pioperties in 
heterogeneous media is that developed by Ben-Amoz (1976), in which the 
concept of heat transfer, as a whole, is redeveloped starting with an analysis 
of “micro- (or local) conduction” within a “micro-medium”. The collection 
of all these micro-media makes up the entire medium. The macroscopic heat 
conduction equations, and the pertinent boundary conditions, are obtained 
using vanation formalism. This approach was extended Baker-Jarvis and 
Inguva (1985) using green function technique. 

Argyns (1994) desenbed Analytical and numencal models for 
mechanical and thermal properties of plasma-sprayed ceramic coatings as a 
function of microstructure. Although these studies were not directed 
specifically toward FGMs, they are nevertheless relevant from the 
standpoint that plasma-spray is one method for ceramic fabrication of 
FGMs. 
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2.3 Models for Functionally Graded Materials. 

In this section, we review approaches that have been used for 
modeling FGM design, processing, and performance Since processing and 
performance of such materials are generally integrally connected, features of 
both are incorporated into various models We look first at models foe the 
spatial vanation of composition, then at system-type models for design and 
performance, and finally at models developed for a variety of aspects of 
FGM behavior 

2.3.1 Models for spatial Variation of composition. 

FGM models generally require an assumption regarding the spatial 
distribution of their constituent phases Consider, for example, an FGM that 
has two constituents, which we denote as 1 and 2. Assume the geometry is 
one-dimensional, with the x-direction being the direction of the 
microstructural gradient We first treat the local volume of phase 1, c, as a 
continuous function (the volume fraction of phase 2 being 1-c if the 
material is fully dense). These functions can be quite simple, the ability to 
exhibit curvature, both “concave upward” and “concave downward”, being 
desirable. 

Wakashima et al (1990), used the following expression for c. 


c = 



(2 4) 


Where x, and X 2 border regions of pure phase 1 and phase 2, respectively, 
and «is a vanable parameter, the magnitude of which determines the 
curvature of c. The curvature can be made concave upward and concave 
downward, to a greater or lesser degree, by proper choice of n . 


Markworth and Saunders (1995) used a quadratic function 1 e. 

c = ao +a,x + f/2x^ (2.5) 

Where the «;(/ = 0,1,2) aic variable parameters whose values are 
determined by imposing constraints and by optimization process itself. 


/ 
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These simple functions still offer considerable flexibility of choice relative 
to the “shape” of the c function 

Markworth et al (1995) plotted the equation 2 4 for various values of 
/land eq 2 5 for the case c(o) = Oand c{x) = l, where X is the thickness of 
FGM along the graded direction These two condition were used to eliminate 
floand < 3 , from equation 2 5, leaving a^as the sole vanable parameter. They 
observed that equation 2 4 offers a greater range of composition profiles 
than does equation 2 5, although the latter equation can be shown to be 
simpler to deal with, from an analytical point of view, than the former. 

2.3.2 Systems approach to FGM modeling. 

A useful approach to the overall modeling of FGM processing is 
based on a system-analysis strategy that has been employed by vanous 
investigators. Hirano et al (1988), have developed an “inverse design 
procedure” which is illustrated by the flow chart shown m Fig 2.3 Here, the 
structure and boundary conditions are specified initially. Then, several 
combinations of matenals are assumed along with different assumptions for 
the spatially dependent mixture ratio. The temperature and thermal stress 
distributions are calculated for these various combinations, and the 
calculations are repeated until “optimum” conditions are obtained Attention 
IS paid to the use of appropnate microstmcture-dependent matenal-property 
models. 

Another systems-type appioach to the optimization process was used by 
Tanaka et al (1993), and consists of the following sequence of steps. 

1 . Choose an initial composition profile. 

2 Carry out a preliminary analysis of non-stationary heat conduction and 
thermal stress. 

3. At each time step, examine the design criteria 

4. If the design criteria are violated, calculate a quantity known as “ thermal 
stress sensitivity increment”. 

5. Find the optimum composition profile that satisfied the design cnteria. 

6. Repeat the analysis of non-stationary heat conduction and thermal stress 
with the composition profile determined 

7 If the design entena are violated at another time step, return to step 4. 

These systems-type approaches to FGM design are just ordered 
sequences of steps, carried out to ensure that the resultant material will 
perform adequately m its intended application However, it is important to 
note that the term “optimum conditions” used above may not necessanly 
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have a unique definition, but may instead involve a set of compromises 
among various quantities (e g material types, heat flow rate, FGM thickness, 
maximum thermal stress, etc ) m order that the desired application is 
achieved over an acceptably long in-service lifetime 


2.3.3 Studies of Thermomechanical Response of FGM behavior. 

Next, we turn our attention to some modeling studies that were 
directed toward a variety of aspects of FGM behavior. These may vary in 
degree of sophistication, but they shaie, together with the studies already 
described, the common goal of maintaining both the performance 
characteristics of the FGM and its structural integrity 

An early version of an FGM was studied by Stewart et al (1986) for 
intended use as a thermal protection system for spacecraft re-entry. The 
anticipated need was high-temperature capability at one end and good 
mechanical behavior at the other end Two different matenals were bonded 
together to serve this purpose Computation of thermal response agreed well 
with expenment 

Markworth and Saunders used their model (1995) to maximize or 
minimize, as desired, the heat flux traversing the graded direction In so 
doing, they found some unusual behavior, as follows they assumed that 
thermal conductivity, modulus of elasticity and thermal expansion 
coefficient were all given by Voigt-typc estimates, for a FGM that was pure 
ceramic at one end and pure metal at the other end (x=0 and x=X, 
respectively). The ceramic end was taken to be at a higher temperature than 
the metal end. Then for certain values of thermal conductivity, modulus of 
elasticity and thermal expansion coefficient for the pure metal and pure 
ceramic, and for conditions corresponding to the extreme cases of maximum 
and minimum rates of heat flow across the graded direction, the thermal 
stress distnbutions are as shown m Fig 2.4. Clearly, the stress distribution 
for maximum heat flow has a maximum value that lies inside the FGM, a 
fact that would have to be accounted for in material design. Another result 
(one that was expected) was that the optimum design for minimum heat flux 
across the FGM contained minimal metal in the structure, subject to imposed 
design constraints, and vice versa for maximum heat flux. 

Result of finite-element calculations of axial, radial, circumferential, 
and shear stresses m disc-shaped FGMs, were reported by Watanabe and 
Kawasaki (1990). They pointed out two considerations that are relative to 
such calculations. One is to minimize the thermal stress caused by cooling 
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from the sintering tempeiatuie. 'I'he other is to minimize the stress generated 
dunng performance, e g as a thermal bamer. 

Another numerical study of thermal stresses in FGMs, associated with 
both processing and performance, was conducted by Matsuzaki et al. (1991). 
They concluded that one of the important factors affecting these stresses is 
the compositional distribution In addition, they found that the key issue for 
this particular FGM (M 0 S 12 -S 1 C/T 1 AI) is the high thermal stress generated 
during processing, a problem that requires furthers attention 

A finite element model, developed by Miller et al (1993) was 
intended to be a guide m both the design and fabncation of a N 1 AI-AI 2 O 3 
FGM The model was used to estimate residual stress as a function of 
structure of the material It also included effects of thermal cycling. 
Elastoplastic analyses were conducted by Giannakopoulos et al. (1994) of 
cyclic thermal response in multi-layered materials These layers consisted of 
a metal, a ceramic, and an FGM between the metal and the ceramic. 
Analytical and finite element methods were both used It was found that the 
smooth gradation of composition between metal and ceramic resulted, under 
thermal cycling conditions, in a reduction of (a) thermal residual stress, (b) 
accumulation of plastic strain, and (c) stress concentrations at free edges. 
Parametric information was also obtained relative to effects of the FGM or 
thermal response of the layered material occurring dunng processing 

Kawasaki and Watanabe (1987) earned out a finite element analysis 
of thermal stresses in metal/ceramic FGMs with a 1400°C-temperature 
difference assumed to exist across the material. A linear dependence of 
matenal properties on composition was assumed. The vanables studied were 
the interlayer width, fineness of the compositional gradation, and the 
compositional profile. They found, for example, that the maximum axial 
thermal stress decreased strongly with increasing interlayer width They also 
calculated this stress as a function of n [see eq 2.4] and found that a value of 
this parameter exists where the stress is a minimum 

Matsuzaki et al. (1992) earned out an analytical design of a ZrOi/TiAl 
FGM intended to insulate the cooling structure of Scramjet engines. Use of 
an FGM having optimized thermal resistance parameters was found to 
reduce the required mass of coolant to just 35% of that needed for a cooling 
structure without an FGM. 

The thermal shock response to heating abruptly an isotropic 
metal/ceramic FGM was modeled by Zhang et al. ('1994), who found that 
thcrmomechanical coupling exerted a strong effect on this response. One can 
thus conclude that the neglect of such coupling effects under conditions of 
non-uniform exposure and sudden, intense heating needs to be questioned. 
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Transient thermal stresses in an FGM plate were studied by Teraki et 
al. using an elastic-plastic formalism. A micro-mechanical approach was 
used, based on the assumption that local microstructures can be described in 
terms of spherical ceramic (metal) particles embedded within a metal 
(ceramic) host phase at volume fractions that vary m a quasi-contmuous 
manner Transient thermal stresses were calculated for a zirconia/stamless 
steel FGM plate, subjected to cycling thermal load, for both elastic and 
elastoplastic conditions A possible thermal ratcheting effect was observed. 

Jin and Noda (1994) who modeled a crack m a semi-mfimte plate of a 
FGM carried out a study of transient thermal stress intensity factors This 
involves a reduction of the thermal and mechanical problems to two systems 
of singular integral equations It was found that substantial lowenng of stress 
intensity factors could be obtained by appropnate selection of 
thermophysical properties 

A non-destructive method for detecting and evaluating the distribution 
of elastic parameters along the graded direction of an FGM has been 
developed by Kuwahara et al (1992) The method that was used was 
“reflective impulse response”, which includes determination of the acoustic 
impedance profile (Acoustic impedance is the product of density times 
speed of sound) The composition profile that was assumed is that given by 
eq 2 4. 

Levit et al (1996) determined the residual stresses in Plasma-sprayed 
ZrOi-lYiOi thermal bamer coatings deposited on Ni-base superalloy 
substrates which were kept at different temperatures. The stress 
measurements were carried out on as-sprayed specimens with X-ray 
diffraction equipment using the “sin^'^” technique. It was found that the 

residual stresses in ceramic layer changes from tensile to compressive with 
increasing substrate temperature during deposition from 75 ‘’C to 500 ”C. 
From microstructural observations and residual; stress analysis it followed 
that the high tensile stresses are responsible for the “segmented structure” 
formed within the ceramic layer, while high compressive stresses lead to the 
development of microcracks parallel to the coating-substrate interface. A 
mathematical model was developed for predicting stress distribution in as- 
sprayed coating-substrate composites. The experimental data showed good 
agreement with predicted result 

Erdogan and Ozturk (1995) considered the antiplane elasticity 
problem for a functionally graded coating bonded to a homogeneous half 
space. The coating is assumed to contain pci iodic cracks pcipcndicular to 
the surface, dlie problem is formulated in terms of an integral equation with 
strongly singular kernels Three dimensionless parameters representing the 
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crack depth, the matenal nonhomogeneity and the crack periodicity are 
identified In addition to the mode 3 stress intensity factor calculated by 
varying these three parameters, the results presented include a qualitative 
discussion of the question of fracture instability, the effect of penodic 
cracking on the relaxation of stresses on the coating surface, and the 
companson of the total strain energy released as a result of surface cracking 
with that assumed in a simple stress relaxation model 

Ostoja-Starzewski et al (1996) considered a unidirectional flber- 
remforced composite with an interphase between the fiber and matrix taken 
as a graded zone of two randomly interpenetrating phases of these materials. 
In particular, they take this mterphase as a functionally graded matenal The 
objective was to present a micromechamcs based method to treat FGM and 
to calculate the effective macroscopic properties (transverse conductivity, or, 
equivalently, axial shear modulus) of such a composite matenal. This 
problem requires the treatment of several length scales the fine mterphase 
microstructure, its mesocontinuum representation, the fiber size, and the 
macroscale level (of many fibers) at which the effective properties are 
defined. It is shown through an example that a convergent hierarchy of 
bounds on the effective response is obtained with systematically mesoscale 
resolutions. 

Finot ct al. (1996) studied experimentally and numcricallly the 
elastoplastic deformation chaiactcristics of a plasma-sprayed, tn-layered 
composite plate subjected to thermal cycling from 20*^0 upto 800°C. The tri- 
layercd solid comprised polyciysrtallinc Ni and Al 20 i oulct layers and a 2.2 
mm thick compositionally graded N 1 -AI 2 O 3 composite interlayer (FGM) 
wherein the composition vaned approximately linearly along the layer 
thickness The experiment involved m situ and ex situ measurements, 
employing a scanning laser technique, of the changes m the overall 
curvature of the unconstrained plate ansmg from the thermal mismatch 
between the constituent phases. The vanations of curvature, accumulated 
plastic strains and thermal stresses at different locations m the layered solid 
were also assessed numencally with the aid of available continuum 
formulations, and the numencal predictions were compared with 
experiments, wherever appropnate. It was shown that when only small 
plastic strains exist in the Ni layer or m the Ni-rich end of the FGM layer, 
known formulations are capable of providing approximate predictions of the 
cyclic vanation in curvature, the onset of plasticity and some features of the 
initiation of cracking (and its location). They also examined possible sources 
of error in the experimental measurements of curvature and in the 
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inteq)retation of thermally induced deformation due to some processing 
conditions 

Taya et al (1996) pioscd a new dislocation punching model for a 
functionally graded material subjected to a temperature change, using 
Eshelby’s model. FGM, consisting of several layer, is deposited on 
aceramic substrate Two types of microstructure were examined for a layer: 
one consists of a metal matrix and ceramic particles and the other of a 
ceramic matnx and metal particles. An elastic energy was evaluated when 
plastic strain, m addition to thermal mismatch strain, is introduced in the 
metal phase. The work dissipated by the plastic deformation was also 
calculated From the condition that the reduction in the elastic energy is 
larger than the work dissipated, a critical thermal mismatch strain to induce 
stress relaxation was determined The magnitude of the plastic strain was 
also determined, when the relaxation occurs. The theory was applied to a 
model FGM consisting of mixtures of Pd and AI2O3 on an AI2O3 substrate. 

Weissenbek et al (1997) analyzed the elasto-plastic deformation due 
to thermal and mechanical loading of layered metal-ceramic composites with 
compositionally graded interface using detailed finite element models. The 
madel matenal considered was a N 1 -AI 2 O 3 layered system with a 
compositionally graded interface. The analyses considered planer geometries 
with perfectly periodic arrangements of the constituent phase, by recourse to 
new classes of square-packing and hexagonal-packing unit cell formulations 
for the graded material. They also considered graded graded phase 
arrangements within which large numbcis of niiciostructural units of the two 
phases are randomly placed within the unit cell. It was found that square- 
packing arrangements provide the best possible bounds for the thermal 
strains and coefficient of thermal expansion of the graded multi-layer, 
among the different unit cell models examined, however, no unique bounds 
could be identified for mechanical loading The numencal predictions of 
thermal and mechanical response were compared with those provided by the 
mean-field approach involving an incremental Mon -Tanaka analysis and by 
the simple rule of mixture approximations. The former method provides a 
stiffer mechanical response than the finite-element unit cell models. The 
finite element prediction of bending coefficient thermal expansion due to 
thermal excursions match the overall trends observed experimentally for the 
N 1 -AI 2 O 3 graded system, and further provide a quantitative prediction of the 
temperature for the onset of plastic flow m the graded matenal. 

Kesler et al. (1997) proposed an experimental method, which enables the 
determination of processing -induced intrinsic stresses, elastic modulus, and 
the coefficients of thermal expansion of surface coatings of homogeneous 
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and graded compositions. In this method, a number of identical substrate 
specimens are coated simultaneously with surface layers of fixed and graded 
compositions, and the specimens with different layer thickness were 
penodically removed from the deposition chamber. It was shown that the 
following result could be obtained from knowledge of the strain or curvature 
and thermal history of the coated specimens, in conjunction with simple 
four-point bend tests and thermal loading/cycling at different temperature: 

1 the magnitude of the processing-induced mtnnsic stresses through the 
thickness of the coatings, 

2 the m-plane Young’s modulus, E, as a function of the coating thickness, 

3. the coefficient of thermal expansion, a, as a function of coating 

thickness, 

4 the vanation of E and a as a function of temperature at any thickness 
location within the coatings, and 

5 the separation of internal stresses ansing firom the deposition process (so 
called “intrinsic” or “quench” stresses) 

The thermo-mechanical analyses underlying this method were discussed in 
detail, and its significance and limitations were addressed The proposed 
method was used to determine the evolution of processing-induced stresses 
dunng the successive build-up of plasma-sprayed N 1 -AI 2 O 3 coatings of 
homogeneous and graded compositions. 

Suresh et al. (1997) presented computational and experimental results 
on the evolution of stresses and deformation fields due to indentation from a 
rigid spherical indentor on a giadcd subshalc. The analyses addressed the 
variation of Young’s modulus, E, of the substrate as a function of depth, z, 
beneath the indented surface for an exponential law, E - E^e°" , where Eq is 

the Young’s modulus at the surface and 1/a is a length parameter. The 
theoretically predicted force-mdenter penetration (P-h) curves were also 
compared with direct expenmental measurements made on compositionally 
graded N 1 -AI 2 O 3 and TiAl-Y 203 -stabilized TZP composites of known 
composition gradients. A new method was proposed for the estimation of 
Young’s modulus vanations through a compositionally graded layer by 
recourse to sphencal indentation 

Dao et al. (1997) developed a physically based computation 
micromechamcs model to study random and discrete microstructures in 
FGMs. The influences of discrete microstructures on residual stress 
distributions at gram size level were examined with respect to material 
gradient and FGM volume percentage ( within a ccramic-FGM-metal three- 
layer structure). Both thermoelastic and thermoplastic deformation were 
considered, and the plastic behavior of metal grains was modeled at the 
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single crystal level using crystal plasticity theory. The results were compared 
with those obtained using a continuous model, which does not consider the 
microstructural randomness and discreteness, yin an averaged sense both the 
micromechamcs model and the continuous model give practically the same 
macroscopic stresses, whereas the discrete micromechamcs model predicts 
fairly high residual stress concentration at gram size level with only a 300°C 
temperature drop in a N 1 -AI 2 O 3 FGM system. Statistical analysis showed 
that the residual stress concentrations are insensitive to material gradient and 
FGM volume percentage The results obtained provide some insights for 
improving the reliability of FGMs against fracture and de-lammation. 


28 




Ipmbinalions 



Optimum Combination and Distribution 
Profile 


•ig 2 1 Mow chail ioi the “invciso design pioccduie” used loi hGM design 
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Fig 2 2 Theimal stress distribution across a metal-ceramic FGM designed for maximum 
heat flow (-) and minimum heat flow ( — ) acioss the graded direction The stiess is 
normalized with respect to the maximum thermal stress in an all metal layer 
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CHAPTER 3 


Temperature Distribution and Heat flow rate. 


The present study is concerned with thermomechanical response of a 
FGM plate subjected to temperature vanation across its two surfaces. In this 
chapter we look at the steady-state temperature distribution across the plate 
as a function of vanous matenal parameters such as thermal conductivity 
and composition profile. In some applications FGM is used as insulator i.e 
to minimize heat flow through structure, m other way one may like to 
maximize the heat flux. The basic equations are outlined in this chapter for 
heat flow rate and temperature distnbution m FGM. Further an expression 
for heat flow rate m bimetallic plate having the same volume of individual 
constiuent as in FGM is also derived. As a specific example, results for a 
given matenal pair (AI 2 O 3 - Ni) are presented and some general observations 
are made. 

3.1 Temperature Distribution in FGM 

In this study a metal-ceramic FGM having a plate like structure is 
assumed. The size of plate is large compared to its thickness. Temperatures 
on ceramic and metal side are and respectively (Fig. 3.1). On the left 
hand side there is pure ceramic and on the right hand side there is pure 
metal The concentration of metal vanes according to a power law from left- 
hand side to nght hand side. The use of power-law expression covers a wide 
range of composition profile by adjusting the value of the exponent n as 
shown m Fig. 3.2. Temperature r, is greater than Temperature T^. and ^2 
are the thermal conductivities of ceramic and metal respectively. The 
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temperature at an arbitrary location x: is given by T at steady state Similarly 
IS the thermal conductivity at an arbitrary location x . 

Following assumptions are made 

1. The y and z dimensions of FGM are large compared to the plate 
thickness that is X 

2 Concentration of metal m FGM vanes only m the x -direction and 
follows power law that is given by the following equation: 

( 

\^ ) 

where n is an arbitrary index and c is the concentration of metal at 

location X. 

3 Thermal conductivity of FGM at an arbitrary location is determined using 
Voigt-type estimate (Rule of mixtures) 

4. The heat flow is essentially one-dimensional in nature. 


Defining 

7] = — , as a non-dimensional distance, we obtain • 


c=n" (3 2) 

where n ranges between 0 and 1 . 

For 1-D steady state heat conduction [Holman, 1996], assuming unit area 
for heat flow: 

where, is the heat flux, T is the temperature, and is the thermal 

conductivity at an location ;c . 

Using rule of mixture for 
k^ =A:,(1-c)+A:2C 

= k,+ik,-k,}f (3.4) 


Defining non-dimensional temperature as 

T -T 
e = i-ii- 

T.-T, 


(3.5) 


Then, G ranges between 0 and 1. 
From equation 3.1 and 3.5, we obtain 

dx X dr] 
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ciT 

Substituting for k , and — - in eq 3 3 

d\ 

Q. = - + {K - -^1 >? " J ^ , which gives 

X dr\ 

^ > 

dr, '(T,-T,)k,(l + mr,-} 

Here 


On Integrating eq. 3 7 from t] = 0 to 1 


jdB =-J 


q,X 




dr] + C| 


e, -e. = - 


oi+"m 


p, , 

f + Qj 


where 6^ is the value of 0 at rj . 

Imposing the boundary conditions 
0=1 at 7j = 0 
0 =0 at T] = 1 


m eq. (3.9) we obtain: 

Cl =0 

and 

- {T,-T,)k, 1 

^ 'f ‘in 

il + mri” 

Substitute these values m eq (3 9) 


(3.7) 

(3.8) 


(3.9) 


(3.10) 


t 
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(3.11) 


9=1 


h 


dr[ 


0 i + 

V Ch) 

Jl + WT]" 


The actual temperatures are therefore obtained by 

/ 

T = T,-(T,-T,) 


r QH 
oi+"”r 




d^ 

1 + WI]" 


(3.12) 


Equation (3.9) gives the value of temperature at any value of tj (o :Siii l). 

3.2 Heat flow rate in FGM and Bimetallic plate, having equal 
volume of metal and ceramic 


(A) Heat flow rate in FGM 


Heat flow rate in FGM is given by eq. (3 10) for a given values of m and n 

„ - (r,-r,K 1 

Jf y dn 
ll + mr]" 


This may be expressed as a non-dimensional parameter H given by 


(r.-r,)*, 


(3.13) 


1 

V dl^ 

•l + mri" 


(3.14) 


(B) Heat flow rate in a bimetallic plate, having equal volume of metal 
and ceramic as in a corresponding FGM. 

Metallic volume fraction m FGM for a given value of n is determined in 
the following manner: 

= Jn'rfn 

0 
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1 

// + 1 

Since 

v,,+v^=\ 

Therefore, 


yc= 


n 


(3.15) 


where is the volume fraction of ceramic in FGM for a given value of n 
Now consider a bimetallic plate having volume fraction of metal and 

ceramic equal and the relative thickness of the layers will be — ^ 

M +1 


and respectively. 

The temperature at the interface of ceramic and metal will have an 
intermediate between 7, and T^. 


The following assumptions are made 

1 . y and z dimensions of plate are large as compared to X. 

2. Area of heat flow is taken unity 

3 Nature of heat flow is one dimension steady state conduction flow. 
For steady state heat flow 



(3.16) 


where is the equivalent thermal conductivity of the bimetallic composite. 

Since heat flux should be same through the metal and ceramic plates (m the 
absence of any heat source or sink). 


^ - K(T,-T,) _ k,{T,-T,) 

« + 1 « + 1 
Therefore we obtain: 


(3.17) 


° nk^ + k^ 


(3.18) 


From eq. (3 8) 
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m == 


^2 ^1 

K 

= {m + 1)^-, 

kj — fk^ 

Where 

r = m + \ 

Substituting k^ = rk^ in eq (3 1 8) 
r, +nrT\ 


T 

■* n 


nr + l 


From eq (3 1 7) and (3 20), we obtain 

t,(„+iX7;-r,> 


7. 


X{nr + l) 


and substituting m eq. (3,16), we obtain 

k 

nr + l 

*„=*.(!+ — i-y) 


« + ! + ■ 


/n 


Substituting this value of k into eq (3.16) 


(3 19) 


(3.20) 


(3.21) 


A,(r,-r,) 


X 


i+- 


n + 1 + ■ 


mj 


(3 22) 


In term of the non-dimensional parameter (//), defined in (3.13) we obtain: 
r ^ 


= 




-T,) 


1 +- 


1 


7J4-1H 

V mj 


1 + - 


1 


n + l + — 

V mj 


(3.23) 


(3.24) 


Note that non-dimensional heat flux has the same definition as for the FGM. 
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Figure 3 1 A metal-ceramic FGM in which concentration of metal varies in x- 
diicction accoiding to powei law (c = /J,") 



Figure 3 2 A metal-ceramic bimetallic plate, having equal volume fraction of metal and 
ceramic as m FGM for an arbitrary value of n 
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3.3 Results and Discussion 


In this section results are presented based on the theory developed in 
preceding sections MATLAB Simuhnk is used to determine the temperature 
distribution and heat flow rate in FGM as well as in the corresponding 
bilayered plate having the same fraction of individual constitute as m FGM. 

Figs. 3.3-3.6 show the variation of non-dimensional temperature 0 in 
the plate thickness direction for various values of composition profiles 


(n= 25-3.0). The value of m 


( 


m = 


k^-k, 




ranges from 1 - 4. 


In each of these cases, there is a smooth change m slope of 0 verses tj 
curve unlike the case of a equivalent bi-layered plate where there is a sudden 
change m slope across the interface This has an important beanng on the 
stress field m the matenal, as described inn chapter 4. 

The vacation of non-dimensional heat flux with composition profile 
and its companson with corresponding heat flux m equivalent layered 
matenal is shown m Figs 3 7-3.12 

For small values of n, FGM is better than the corresponding bi-layered 
matenal in maximizmg/removmg heat. However, there is a cross over point 
(ncr) , beyond which the bi-layered material will be more effective If the 
intended application is thermal insulation, the FGM is better for values of n 
greater than n^ . The position of cross over point is dependent on the ration 
of the thermal conductivity of the two matenals as shown in Figs 3.7-3.12. 

The vacation of ncr with thermal conductivity ratio has been obtained 
through a parametnc study and is shown in Figs 3 13. 

From fabrication point of view, it would generally be difficult to make 
matenals with n>10 or so, therefore a detailed view of Fig. 3.13 in the range 
n=.l- 10 is given m Fig. 3.14. 

This information serves as useful design guideline in selecting the 
constituent materials and specifying the appropriate composition profile to 
maximize the performance of the structure. 
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Fig 3 6 1 Temperature distribution in Bi-layered structure, which corresponds to 
a FGM, having the value of n as 1 (n=l) 
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Varialion of H with niot m»1 5 
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Fig 3 9 Variation of heat flow rate in terms of H (non-dimensional parameter) in 
FGM and bimetallic plate for kj=2 5k^ 
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m FGM and bimetallic plate for k^ = 3/c, . 


45 



Variation of H with n for m«3 



Fig 3 II Variation of heat flow rate m terms of H (non-dimensional parameter) 
in FGM and bimetallic plate for = 4^, 
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Fig 3 12 Vanation of heat flow rate in terms of H (non-dimensional parameter) 
tn FGM and bimetallic plate for = 5^, 
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Fig. 3. 1 3 Variation of ricr with m 


m = — =■ 






for a FGM/bilayered matenal 


47 





CHAPTER 4 


Thermal stresses 


1 his chapter deals w ith thermal stresses in a FGM plate as a function 
ot \arious niatciial paiameters and composition profile when the plate is 
subjected to tcnipciatuie gradients. The temperature distribution required for 
such calculations has already been described in the preceding chapter. 

Attention is paid to determination of thermal residual stresses 
developed dm mg fabiication of FGM. A trilayered ceramic-FGM-metal 
structure is considered & subjected to 1 00 “C temperature drop Such 
structure is often employed for mctal-ccramic joints and is of great practical 
significance. 

4.1 Determination of thermal stresses in FGM 

We consider classic beam and plate theory analysis of stresses and 
curvature, which develop in the layered solid during temperature 
changc[Boley and Weiner, 1960] The resulting kinematics does not resort to 
the beani/platc appioximations but are exact within the context of small 
strain analyses. The model system investigated is ceramic-metal FGM, in 
which composition changes gradually and continuously according to power 
law, from pure ceramic to puie metal For this continuum approach to be 
valid, the thickness of FGM must be significantly larger than the 
characteristics nucrostructural dimensions (e g. gram size). Creep of 
constituent matenals is neglected. The FGM is assumed to be initially stress 
free. 

The geometry of the FGM considered here is such that the problem can be 
reduced to one dimension. Fig. 4.1 schematically illustrates the geometry of 
the FGM. I'he total length, width and height of FGM are 1, w and 2h, 
respectively. When l,w»h, the geometry modeled represents plane stress 
conditions. Here, the variables of interest in the study of thermomechanical 
response depend only on the vertical position x(along the thickness 
direction). 

It IS also assumed that thermomechanical and physical properties follow 
voigt type estimate . Only elastic response is considered 
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Composition of mttal vary according to power law(eq 3 lb), that is 

^ = ( 41 ) 

hei e 

(4 2) 

Define 



(4. 3) 


Now relationship between jj and 
transfomiation relationship. 


V = x + h 
■ 


X _ 

2h~2 


+ 1 




P IS determined using coordinate 


h='^(l + /i) 

Substituting this value in eq (4.1) 



(4.4) 


(4. 5) 


£, and £2 are the elastic modulus of ceramic and metal respectively. 
Similarly a, and are the thermal expansion coefficient of ceramic and 
metal respectively. 

Now, using rule of mixtures (voigt type estimate) for above properties 


E = E^{l-c)+E,c 


E = E,+(E,-E^)c 

(4,6) 

a =a,(l -c) + a 2 C 


a =a, +(a 2 -ai)c 

(4.7) 
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In gencial, the components of any three dimensional strain field arc 
subjected to six conditions of compatibility [Snnath, 1996] But if the strain 
IS allowed to be a tunction of x only, that is in the present study, then four ol 
these conditions aie satisfied identically and the remaining two condition 
both reduce to 


e(0 = Cn + (4 6) 

Here «is the strain atx = 0 and b is the curvature of the beam in its 
plane. In plane sticss condition, the only non-zero stress component is 

a(x) = cr„. 

Therefore 

c7(x)=£(x|f.(x)-a(x)Ar(x)] 
a(x)= £(x)[a + Zjx -a(x)Ar(x)] 

cr(.x) = «£(x)+/j.x£{.x)-£(x)a(x)Ar(x) (4. 7) 

AT represents the change in temperature from initial stress free state 


1'he resultant force and the resultant moment of the stress distribution 
a(x) along the height, x, must be equal to zero as there is no external applied 

force and moment. 

Using force balance: 

h 

jCT(.t)<ilx = 0 

-A 

Put the value of cr(x) from eq (4 9) 

h h h 

a j£(jx:)iic + b |x£(x)fix = j£(x)a(x)Ar(x)ix: 

-A -A -A 

oAj + bA^ — Aj 

Using moment balance: 

A 

|a(x)xxix: = 0 

-A 

Substitute the value of ct(z) from eq (4, 9) 

AAA 

a ^zE{z)(k + b jz^E{z){Iz = |z£(z)a(z)Ar(z)(/z 

h H -h 


aAj + bA^ 




(4.9) 
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On soh mg eq (4. 1 0) and (4 II), we get 

"'oi7-4,‘' 


j,-’- ot. 

Heic 

j, . 

-h 

h 

J, = |\£(i)/\ 

h 

h 

|£(\ Jti(x)A J(\)£Zx 

-A 


A, = |x'£{x)£x 

-h 

h 

A , = Jx/:(x>4x)A7(x>/x 

-A 



(4 12) 
(4 13) 



(!x=h(l(i (4. 14) 

Limits of p varies between -1 to 1, corresponding to limits of x which varies 
between - h to h . 

Put these limits and c/x = help in above expressions of A^,A^,A^,A^and A ^ . 


1 


if 

(4. 15) 

= h^]pE{p)llS 

(4. 16) 

= h]pE{p)!x{P)AT{l3)ip 

-1 

(4. 17) 

=e]p^Eip)cip 

t 

(4. 18) 

^iPpEiPHfiy^HPW 

(4. 19) 


„I 
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4.2 Condition tor zero thermal stress on ceramic side 

1 hermal stress iii I'CIM is given eq 4 9 which is as follow 
(t(x)= a/;(x) t /n/;(\)- £(\)tx(x)Ar(x) 

01 

cy(/l )= ) f />/!/■ (/n- E{li)a{PhT{(i) 

For ceiannc suilacc, 

= and Ar = (ArL, 

Subslitutiiii^ tile St \ ciliits in sbovc ccjuHtion 2,nd ccjuEting it to zero, wc get 

(4 20) 

This expicssion is independent of a, and because a and 6 contain 

these term, hence these stands canceled out in the above eq 4.20 From the 
above condition we can find out the value of a and b, which will give zero 

stress on ceramic surface The value of a and 6, depends on and 

a, ’ Aj 

E (X li 

n . So for given values find the composition profile n for 

i S I (J! J A J 

which Ihcie is zero stress on ceramic surface. 


4.3 Condition for zero thermal stress on metal surface in FGM 

Similarly a condition has to be satisfied for having zero thermal stress on 
metal surface, which is as follow 

«+/> = 0 (4.21) 

Remaining argument remain same as in section 4.2 

4.4 Determination of thermal residual stresses in FGM 

The model system investigated here is a three-layered composite, where 
the two outer layers comprise a metal and ceramic; the in-between interlayer 
is compositionally graded wherein the composition changes gradually (and 
continuously) from pure ceramic to pure metal. In an attempt, to illustrate 
quantitatively the salient features of the icsull and to provide specific 
examples, we consider Ni and AI2O3 as the model metal and ceramic. 
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rcspccti\ cly Ihc I‘GM layer is a composite of the Ni- AI 2 O 3 mixture, with 
varying piopoitions of the two phases. Same assumptions regarding the size 
of the plate and material piopcrties have been made here as in previous 
section 4 1 

I he I'tiM layei extends fiom x=-a to x = +a It is assumed in this study 
that a = 6 h Only elastic response is considered 

Formulation of the problem is same as in section 4 1 The only difference is 
this that, heic, \r = --lOO‘’C(diop) 

Heic composition of metal in FGM is assumed to vary as follow: 

tor -a < X <+a 



Define 


-1 = ;. 

a 


(4. 27) 
(4. 28) 


Then 
(h = adk 
Now 



Using rule of mixtuics foi E and a . 
E = E,{\~c)^E^c 


a =a,(l -r) + a^c 



(4. 29) 

(4. 30) 

(4.31) 

(4.32) 


From section 4. 1 

h 

-A 

-an h 

A, = |/v/jr-i j£(x>A 1 J/ijc/t 

h » 
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(4. 33) 


(4. 34) 


(4. 35) 
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f, -/'.i ur, 

' \ 3 1 

SubsUlLiie 
X = aX und i/,x " (ill? 


1 f W X 

~3 j"^ 


1 x'f,w ^ 

‘ ' + a’jA-£(A>a 


A, = Jxa{x}/:(x)A7(xVx 

h 

Heie 

Ar = - li ) o ' c ' 

r r "r * 

-J, =-10()j Jxa,/;,f/x+ Jxa(x)£(x)+ Jvaj^jcic 

L ^ fl a 

Substitute 
r = aX and dx = «(//. 

4 = _100j ( iL;£ll(a,£, -«,£,)+ A' ]xa{x)E{X)dX 

-1 


(4. 36) 


(4. 37) 


From section 4. 1 

cr(x)= a£(v) + />\£(x)- «(x)£(x)Ar(v) 

1 Icrc 

Ar = -ioo"c 

Hence 

a{\) = aE{x) >t bE{\)ih \ li)i)a{x)E{x) 

The values of £{a) and a(x) are different in different layer. So the 
expression of ct(x) will be different m different layer. The expressions are as 
follow: 


cr(x) = a£, +/h£, +I00a,£, 

for 

-h^x<-a 

(4. 38) 

a(v) = «£{x) + /n£(x)+ 100a(-x)£(v) 

for 

-a<x£a 

(4. 39) 

cr(v) = a£2 +/>x£2 -1 lU0aj£2 

for 

a<x^h 

(4 40) 
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4.5 -l liermal stresses in bimetallic plate, having equal volume 
of metal and ceramic as in FGM 

For compa, ,son of metal-coam.c FGM with bimetallic plate, having equal 
\olunic liaclioiuil metal and ceramic, it become necessary to determine the 
thernial stic.s.s dislnbution in bimetallic plate using same assumption and 
conditions as vcie m I'-GM thermal stress study. For determination of 
thcnnal .sticsscs in bimetallic plate, it is necessary to know the temperature 
distnbulion in bimetallic plate, 

(A) Temperature distribution in bimetallic plate 

Fiom section 3.2, \olume fiaction of metal and ceramic in FGM is given by 
the tollouing expicssions 

= — i- 

" /Hi 

and 


‘ « f I 

Now consider a bimetallic plate, which has volume fraction of metal and 
ceramic equals to and i;, icspcctivcly 

From section 3.2 
^ 7 ; 1 n{m 1 l)/’j 
*' M(/fMl)ll 

Hence tempeialuic at an aibitrary location tj in ceramic part is given by the 
following cxpicssion 

r, = r, -u-p o<n<~ (4.34) 

/i + l 

Where u- is a constant. 

At tj = -1 - , 7 ; = 7„ 

n f 1 

Putting this condition in equation (4. 34), we get 
('■->■1X7', -r„) 

n 

Put this value m eq. (4. 34) 

(4.35) 

n 

Similarly the temperature at an arbitrary location fi'in metallic part is given 
by the following expression: 
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j.=n- £1 

?r + 1 

Where > is a constant 

A.t t] -- ' ■ , - /|) 

n r 1 

Putting this condition in cq (4 36), we get 

y = ih~7])i,nt} 

Put this \ aluc of 1 in eq (4 36) 

7; = A 


(B) Thermal stresses in bi-Iayered plate 


From section 4 2, 

I] = -^ {1 + /j ) 

Wheie 

j] =-- and /i 
2 b h 

Now 

Tj =o=> /? =-l, 

n .. n- \ _ 
4 = — “=>/^ ; 

;i + 1 /Hi 

7] = I => fi ~ ; 

1 lencc 


2/7 


T = T + 


(/, fiX7;,-7;Xi-;i) 


Now from section 4. 1 
/(, =A jM/tX//! 

/" « 1 

I 

I + j£,4/3 


-I 




2b 
n + 1 




-I 


(4. 36) 


(4. 37) 


(4. 38) 
(4. 39) 


(4. 40) 
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! n I 
«* I 


= /r 


I ■. J/i/-.,//} 


-I I 


:nA 
(n 4 1 ) 




I 

n * 1 t 

\a,E,{\T\dp^ la,E,{^T)„dli 


/L = h 


-I 


« j 

»t + ! 


Here 

(Ar), = r, - 7 ; 

Wheie r, IS assumed as room temperature. 
(a 4 




Similarly 

(ATI = T'. - n 

(^r) JnA]U,-'Ti\-li) 

Hence 


= // 


(4.41) 


/4 ' j I ( 7 ; _ r, )- (;i±%ZMU j 


I 


2n 


n-\ 

'nTl 


(4 42) 


1 

"I ji^E^dfi+ jp^Ejdji 


A. =h 


n I 
«+! 


= 3 Tu" [^ 1 1'* - 0’ + (« + iF )+ ^2 I" ■'■ ~ ~ 0^ 1 

3 (/i + 1 j 


(4. 43) 


/l, =//^ J/iAt/})u(/l )A7t/J>//i 

-i 
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( 4 . 44 ) 


Expiessions tor thermal stress will remain same as in 
\ allies ol a, b, and expression for AT will change. 


previous section. Only 
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4.5 Results and Discussion 

Model 5 .> stem elioseti for this study is a N 1 -AI 2 O 3 FGM The choice of 
tins model s> stem ,s pt edicaled upon the following considerations' 

1 . Ilie tlieimonieehunical and physical properties of Ni and A1,0, are well 
known as these materials have been widely studied 
2 Compositionally graded microstructures based on Ni and AIoOt have 
been pioduced using a vaiiety of available techniques. 

3. Ni-bascd alloys aie extensively used is aiicrall turbine engines where the 
thcimal cycling response of the alloy bonded to ceramic coatings is a 
topic oi considciable intcicst. 

4. Ni, has a high inciting point (1726 K), thcic exist many practically 
significant tcmpciatuie regimes, well below the melting point, wherein 
the thcrmomcchanical response can be studied without incorporating the 
complicating effect of creep, diffusion and chemical reactions 

I'he numerical values of material parameters used are as follow: 

E, = 380 CJpa 

£, = 214 ( jpa 

a, =74x10 "'A ' 

=15 4<10 '"‘A' ' 

Thermal stress aic detcrinined m FGM for r, = 120 °C and = 20 °C(room 
tempeiaturc) by using MATI.AB simulmk for different composition profile 

and " j • 

Figs. 4.3, 1 - 4.3. 1 i show the variation of stress along the plate thickness as a 
function of n where n ranges between 0 75 to 4 and m=\. 

The main points to note are: 

(a) The stress variation is continuous, unlike that for a bi-layered material 
where a discontinuity exists at the interface. Sudden changes in stress 
across a layer are potentially damaging for the matenal and must be 
avoided. 

(b) The stresses on the ceramic side may be tensile or compressive 
depending on the value of n . Since ceramics are bnttle in nature and fail 
through crack propagation, it is desirable to have compressive stresses, 
which inhibit fracture. Thus, the present analysis offers a valuable 
approach to enhance the reliability of a FGM through proper choice of 
composition profile (value of n ). 

Similarly results for m = 2 also have been obtained (Figs, 4 4,1- 4.4.1 1 ). 
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The MiiKilion ett thermal residual stress with position for a tn-layered 
mateiial uheie the iiUei mediate layer is a PGM and the surface layers 
conespoiid to puic ceiamic and puie metal respectively are presented m 
logs 4 5 1- 4 5 ^ lliese lesults aie obtained for a drop in temperature of 

100 (' 

The mam obsci\ alums aie as follow. 

(a) 'I’hcic IS no shaip change m stiess value at the two interfaces, although 
the siicss \ anation is not continuous for n = l . 

(b) The discontinuily in the slicss vaiiation at the interfaces can be removed 
thiough piopci choice of composition profile. For example, for n=Q. and 
3, the slicss changes continuously at the ccramic/FGM interface. 

Smulaily, the slope at the second inleifacc may be adjusted. Although the 
FGNh metal inteilacc is mure forgiving 

Depending on the intended application, it may be desirable to have zero 
stress on the cciamic or metal surface through proper choice of composition 
profile, 'fhis is shown in Fig. 4.6- 4.9. 

In view of the buttle natuie of ceramics, it is desirable to minimize the 
tensile stichses on tlie ceramic side Accordingly a parametric study is 
earned out to detcimme the composition profile («*) which would give zero 

( E ^ ( ct k 

stresses on ceiamic side as a function of 5 5 =-^ , y r =-7 ^nd 
These results aie piesented in Figs. 4.10-4.22. 
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H* 




Fig 4.2; A schematic representation of ccramic-FGM-metal 
plane stress model (a=0.6h) 
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Fig 4 2 I Schcnuitic illustration of a bi-Iayered material, which have the same volume 
fraction of indii. idual constituents as in FGM having the composition profile n 



F,g4.2 2: Vanalion of thermal stress rn b.-layered structure, whroh ts equrvalent to FGM 
having the value of n=l . 
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0 06 1 


beta > 

ation of thcmiai stress m FGM for n=0 75 and m=l 


^4? I* cr» thermii sirtsj in FGM for n«1 andm-l 


Fik 4.3 6' Prediction of variation of thermal stress in FGM for n-1 50 and m-1 



ikmrmM • 



beta > 


Fig 4 3 8 Prediction of vanation of thermal stress in FGM for n-2 and m 1 . 
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6 


Variation of thorma! stress in FGM for n-2 6 and m«1 


T 


n— 2.6 



ediction of variation of thermal stress m FGM for n=2 50 and 




Ihenmai str«ss fWpa) 


Variation of thermai stress in FGM for n»»4 and m«»1 



Fig 4 3 11 Prediction of vanation of thermal stress in FGM for n=4 and m- 


















Fig 4 5 2 Prediction of the thermal residual stress distribution in ceramic-FGM- 
mctal tnlayered material for AF = -lOO^’C and n=2 



Fig 4.5.3 Prcdiclion ofllic Uicnlial residual stress dislribuUoii m ccramic-FGM- 
metal trilayered matenal for Ar = -100 C andn-3 
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Fig 4 6 Illustration of zero thennal stress on ceramic surface for m=l Corresponding 
value of n is 1 628 



Fig. 4.7* Illustration of zero thermal stress on ceramic surface for m 2 
Corresponding value of n is 1 21 
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Fig 4 8 Illustration of zero thermal stress on metal surface for m=l The 
corresponding value of n is 3 5166 



Fig. 4 9 Illustration of zero thermal stress on metal surface for m-2. The value of 
nisL7516. 






Fig 4 1 1 Surface view of variation of n* as a function of y 
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Fig 4 12 Surface view of variation of n* as a function of y 
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CHAPTER 5 


Effect of Microstructure 


5.1 Thermophysical properties of heterogeneous materials 


As already has been noted, effective Thermophysical properties of 
heterogeneous matenals are intimately related to the microstructure of the 
matenal Along these lines Nan (1993) has identified three basic type of 
geometne morphologies associated with the distribution of phases within a 
two-phase material 


( 1 ) Dispersed grain structure This is characteristic of a matenal for which 
the volume fraction of one phase is low, and is discretely and randomly 
distributed within the host phase [Fig. 5 la]. It is also characteristic of a 
material for which the volume fraction of the dispersed phase is 
somewhat higher, but the phase remains discretely and uniformly 
dispersed [Fig. S.lbj. 

(2) Aggregated gram stiuciure. This morphology characterizes the case for 
which the volume fraction of the minor phase is increased to the point 
that it is no longer discrete, but is aggregated in the form of clusters of 


definite size [Fig. 5.1c]. 

(3) Percolation-like cluster structure As the volume fraction of the minor 
phase IS increased still further, it reaches a cntical value, called the 
percolation threshold, at which it is interconnected to such an extent that 
there exists a continuous random cluster, as well as smaller aggregated 

clusters I^Fig- 5# Id]. 

A great deal of effort has been devoted, over the years, to the development 
of thermophysical-property models for a dilute, discrete dispersion of one 
phase within the other. However, the case of intertwined networks of two 
phases as illustrated in Fig. 5.1d present particularly interesting behavior. If, 
for examples, the electneal conductivity of one phase is significantly higher 
than that of other, then a sharp increase in the effective thermal conductivity 



of t le mixture would take place at the composition for which the relative 
amount of the higher conductivity component reaches its percolation 
threshold value A schematic illustration of this behavior is given in Fig.5.2. 
Indeed, a vantty of topological features of the microstructure need to be 
consideied m order to develop physically realistic property models These 
include directionality, connectivity, and irregularity of the phases that are 
present [Markworth et al , 1995] 

A question that must also be addressed is that of how much information 
regarding microstructure must be included in order to describe adequately 
material properties and behavior. Becker and Richmond treated this 
particular issue, as it relates to effects of grains, pores, and second-phase 
particles in numerical simulations of the mechanical response of ductile 
materials. I'hcir models were limited to two dimensions, although the need 
for inclusion of micro-structural features, in order to be physically realistic, 
was demonstrated 

The description of these complex microstmctures can involve some unusual 
difficulties, one example being the case of percolation. Chemikov and 
Rogalsky (1994) have pointed out that studies in two dimensional 
percolation theories, have, over the past decade, borne considerable fhiit. 
However, the three-dimensional problem presents greater difficulties, there 
still being no exact analytical treatments for lattice percolation and random 
continuum percolation, although numencal studies have been earned out. 


5.2 Percolation Theory 

Percolation theory has been studied for various physical (transport) 
properties by a number of workers over the last two decades and 
thecurrent approaches are desenbed in vanous texts such as Efros 
(1986) and Staufer and Aharony (1992). 

The phenomena best described by percolation theory are critical 
phenomena. They are characterized by a critical point at which some of 
the properties of the system undergo abrupt changes. Critical 
phenomena include second-order phase transitions (e.g. the transition of 
a metal from its normal to its superconducting phase when its 
temperature is lowered). The physics of all cntical phenomena si very 
unusual, but there are some common features, the most important of 
which is that in the neighborhood of the cntical point, the system 
appears to break into blockks which differ in their properties, with the 
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size ol the individual blocks growing until the system approaches the 
critical point The blocks are quite randomly shaped In some 
phenomena the whole configuration changes chaotically because of 
thcimal motion, while in other phenomena the configuration may be 
lio/cn in time but changes from specimen to specimen. The blocks are 
m complete disorder, so that no regularity is discernible from an 
instantaneous photogiaph However, this geometry, which can be called 
the geometiy oi disoider, has quite definite properties “on the average”. 
Acluall} , geometiy is inseparable from physical properties For instance, 
the ph) sical piopcities of a crystal are determined by the geometry of its 
lattice Likewise, the “geometry of disorder” determine a number of 
piopeities of a system in the vicinity of a cntical point The most 
interesting feature is that owing to the large size of the blocks the 
geometry' is virtually independent of the atomic structure of the matenal 
and thus possess properties common to a number of quite dissimilar 
systems; hence, the universality of the physical properties that we find 
m the neighborhood of critical points This type of relation between 
physics and geometry can be traced in percolation theory 

5.3 Percolation in heterogeneous materials 

As the volume fi action of the minor phase is increased m heterogeneous 
materials, it reaches a critical value, called the percolation threshold, at 
which It IS interconnected to such an extent that there exists a 
continuous random cluster, as well as smaller aggregated clusters. 

As soon as the volume fraction of the minor phase reaches its 
percolation threshold, there is a steep change in transport properties of 
heterogeneous mixture This change is proportional to (v-vj'’, where 
is the percolation threshold and pls^ constant (Staufer, 1992). 

If thermal conductivity of one phase is significantly higher than that of 
the other, then a sharp increase in the effective thermal conductivity of 
the mixture would take place at the composition for which the relative 
amount of the higher-conductivity component reaches its percolation 
threshold value. This can be demonstrated through a simple numerical 
exercise whcie a two phase microstrucluic is created numerically. A 
grid consisting of specified number of cell is chosen, Ihe number of 
cells corresponding to a given volume fraction is determined ^d 
specific cells corresponding to the second phase aie chosen using 
random number generation. The range of random numbers is chosen so 
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as t*.) co\ei the total numbers of the cells The results of such an exercise 
ai e desci ibcd in section 5 5 

5.4 Percolation in FGM 


I 


I he iinpuitancc of pet eolation in heterogeneous materials has been 
discussed in the picccding section. It had implicitly been assumed that 
the heteiugeneity IS spread uniformly through the material 
In the case ol a I'CIM vve arc dealing with a material, which is 
heleiogeneous yet thcie, is a bias in the distribution of the second phase. 
Percolation in such materials may be analyzed m a similar manner (as in 
section 5.3) except that the effect of bias (concentration profile) must be 
accounted for. A numerical study of the same was earned out as follow: 

A gild conipiising of 45*20 cells was generated The continuous 
c'ariation ot concentiation was assumed equivalent to step-like vanation 
thiOLigh a la}er {20*5 cells each). The volume fraction for each layer 


was detci mined so as to give equivalent linear vaiiation over the entire 
thickness The coiresponding number of cells representing the 2"^* phase 
was identified in each layer and their sequence determined through a 
random number routine. The results of such an exercise are shown m 
section 5.5. 

In terms of incorporating such effect in calculation of temperature 
field and thermal stress, Surersh et al., have considered 
micromcchanical appioach whore they simulate such microstructuie and 
model heat transfer and stress m each gram We suggest an alternative 
approach where the overall microstmcture effect of formation of 
percolation cluster on the transport properties is considered in terms of 
vanation of the said property with concentration (and therefore, spatial 
position m FGM) and these expressions are directly used m conjunction 
with the equations denved earlier (chapter 3 and 4). The vanation of a 
property with volume fraction (incorporating the percolation effect) is 
adequately described by an expression having the following form. For 
thermal conductivity for example: 


Define a nondimensional parameter k as follow 


(4.22) 
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(4.23) 



then 


A - 1 i- (4 24) 

w hci e the nibols used has following meaning 
A, - Mfectne theimal conductivity of FGM at an arbitrary location. 
A, - 'I hennal conductivity of ceramic in ceramic-metal FGM 
A, 1 hcmial conducts ity of metal in ceramic-metal FGM. 

V = Metal \ olume fraction at an arbitrary location in FGM. 
v^= Percolation threshold of metal volume fraction. 

/>= An ai bill ary’ index 

Values of these parameter taken in this study are as follow 
A, --=30 W/m K 
A,=300W/mK 
v*-.20,.25, 40 
5 


Now tempeialuie distribution in FGM is given by the following 
expression; 




<//3, 


30 + 270 


0=1 


\-e 




rf/i, 


30-1 270 


1-c 


(;)' 


(4.25) 


Assuming liner variation of metal volume fraction m FGM, then 
V = i3, . 1 lence above expression become as follow; 
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'*) ^ 270 


0 I 
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1-c " 
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0 ) ♦ 270 


h 


I PA' 


(4 26) 


Siniilaily the expicssion for H (H~ 




il 


1 


i 

1 I 

i I 

1 + ni 



) IS given as follow. 


li.xpi ession for thermal stress distribution remain same as eq. 4.9, with 
change in the values of «, 6 and variation of 0 . 


1 
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5.5 Results and discussion 


! ui 5 ^ depicts niiciostuictures of helciogcncous material for different 
volume 1 1 actions ol individual constituent 

InitialK thcie is no inteitvvincd stiucturc, but as soon as the volume fraction 
of phase 2 is about dO'/o the pei eolation effect appears 
Similaily big 5 4 sluws microstructure of an FGM of two phases 1 and 2. 
Phase 2 is lepiesented by black box It is observed that at about 25-30% 
volume fi action of phase 2, pei eolation effect starts in the FGM, hence a 
steep change in the tianspoit piopcrties of the material is expected. 

Based on equation 4 24, Figs. 5.5- 5 7 shows effect of percolation on thermal 
conductivity \anation in FGM for different values of pci eolation threshold. 
Figs. 5.8- 5.9 shows the conesponding variation of temperature along the 
thickness diiection in the FGM. 

The tenipeiatiuc \anation in FGM with percolation effect (for v, = 25) is 
compaied with conesponding FGM without percolation in Fig. 5.10. The 
mam observation is that with percolation results in sharp decrease in 
temperature m FGM for tj< 25 whereas, the slope of the e/ii curve is almost 

constant for FGM without percolation 

fhe siiinihcance ol this result can be seen in Fig 5 12, where thermal 
stresses aic compaied for a linear variation of composition of metal. It is 
observed that with percolation effect, there is sharp decrease m temperature 
on ceramic-rich icgion, which results in compressive stresses on ceramic- 
side. I'his IS beneficial from design point of view since ceramics are bnttle 


and fails easily under tensile stresses. f ,c 

The effect of the value of percolation threshold on the ceramic surface is 
examined m Fig. 5,11, where it is observed that a lower v. g‘ves h.gh« 
compressive stresses. It is speculated that the percolation threshold may be 
lowered by altering the relative gram sizes of the constituent materials. 
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Fi 5 1 

function of vota™ racnor‘''“" '” “ 




Fig 5.4. Microstructure of a FGM depicting the intertwined network of black box at 25- 
30% volume fraction (black box represent the metal gram) 



nrwtaJlic vo!um« fraction •• > 


I-ig.5 5 Schematic illustration of the vanation composition of the thermal conductivity 
of a two- phases mixture, one being good conductor, the other a poor conductor 
Percolation threshold is .20 



Fib 5 6- Schematic lllustralion of the vanaHon composition of llicthcnnal coiiductivily 

of a two-phase mixture, one phase being a good conductor, the other a poor conductor. 
Percolation threshold is .25. 




metallic volum® fraction — > 


Fig 5 7 Sclicmatic illustration of the vanation composition of the thermal conductivity 
of a two-phase mixture, one phase being a good conductor, the other a poor conductor 
Percolation threshold is 0 40 



Fig 5 8 Prediction oftemperaturedistnbution in FGM with percolation effect 

(percolation threshold is 0 20) 






Fig. 5 10 Companson of temperature distribution in FGM wit and without 
percolation effect Percolation threshold is 0 25 







Fig 5 1 1 Prediction of thermal stress distnbution in FGM with percolation effect 
for various values of percolation threshold 
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Fig 5 12 Companson of thermal stress distnbution in FGM with and without 
percolation effect Corresponding percolation threshold is 0 25 





CHAPTER 6 


Conclusions and Suggestions for Future work 
6.1 Conclusions 

This work has dealt with the thermomechanical response of functionally 
graded materials. A plate like structure subjected to sharp temperature 
gradient in the thickness direction is considered. Theoretical models and 
numencal results for temperature field, heat flux, thermal residual stress, and 
thermal stress have been presented both for functionally graded material and 
equivalent bi-layered materials Significant trends and design guidelines to 
select appropnate composition profile for specific matenal pairs have been 
identified and presented The role of microstructure, especially in regards to 
percolation phenomenon m transport properties and its significance for FGM 
are discussed and a novel approach to incorporate such effects in design 
calculation is presented 

The major findings are' 

1. The steady state temperature field m a FGM subjected to temperature 
gradient along thickness direction vanes in a continuous manner unlike in 
bi-layered material 

2 By proper choice of composition profile, stresses in FGM can be 
optimized to suit specific requirements For example, It is possible to have 
zero stress on ceramic or metal surface in FGM, by adjusting the 
composition profile. 

3. There is no sharp decrease m stress profile across the plate thickness 
unlike as in bi-laycred material. Hence there are lesser chances of 
matenal failure by crack propagation. 

4. The discontinuity m stress profile in ceramic-FGM-metal tn-layered 
structure can be avoided by proper choice of composition profile. 

5. Percolation effect become dominating, at a particular volume fraction and 
at this volume fraction, and steep change m the transport properties is 
expected. 



6 Temperature distribution in FGM with percolation effect is continuous, 
but there is steep change in temperature for v < It is advantageous in a 
way to have moie compressive stresses on ceramic surface m FGM 

6.2 Suggestions for Future work: 

1 Heat flow rate through FGM with percolation effect and without 
percolation effect should be calculated Basic equations have been 
developed in the preceding chapters 

2 Behavior of FGM, in plastic region should be investigated 

3. Nucleation of crack and failure of FGM by crack propagation should be 
investigated 

4. There is a need to investigate Creep behavior of FGM since FGM are 
used in high temperature applications 

5. A thought can be given to FGM having 2-dimensional and 3-dimensional 
composition gradient. 
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